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Abstract 

We discuss the possibility of coexistence of spin density wave (antiferromagnetism) and triplet 
superconductivity as a particular example of a broad class of systems where the interplay of mag- 
netism and superconductivity is important. We focus on the case of quasi-one-dimensional metals, 
where it is known that antiferromagnetism is in close proximity to triplet superconductivity in the 
pressure versus temperature phase diagram. Over a range of pressures, we propose an intermedi- 
ate non-uniform phase consisting of antiferromagnetic and triplet superconducting orders. In the 
coexistence region, we propose a flop transition in the spin density wave order parameter vector, 
which affects the nature of the superconducting state. 

PACS numbers: 74.70.Kn, 74.25.DW 
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The competition or coexistence of magnetic order and superconductivity is a very impor- 
tant problem in condensed matter physics. There is a broad class of systems that present 
magnetic order and superconductivity in close vicinity. One of the most important sys- 
tems are the Copper oxides, where singlet superconductivity is found next to antiferromag- 
netism [lj. Another interesting system is Strontium Ruthenate Sr 2 Ru04, where the prox- 
imity to ferromagnetism has been argued as being important to the existence of possible 
triplet superconductivity in these materials Q]. Furthermore the ferromagnetic supercon- 
ductors ZrZn 2 and UGe2 have stimulated a debate on the coexistence of ferromagnetism and 
triplet or singlet superconductivity j^,^). However, unlike any of these previous examples, 
we discuss in this manuscript a quasi-one-dimensional organic superconductor, the Bech- 
gaard salt (TMTSF) 2 PF 6 , which has a phase diagram of neighboring antiferromagnetism 
(AFM) and triplet superconductivity (TSC) P|. 

The antiferromagnetic state of (TMTSF) 2 PF 6 is present at temperatures T < 12 K and 
pressures P < 6 kbar, and is characterized by a spin density wave (SDW) [5|. The SDW order 
parameter N (Neel vector) is anisotropic, having an easy axis along the crystallographic b' 
axis 6], which is also the intermediate direction for conductivity. This antiferromagnetic 
state is supressed at pressures higher than 6 kbar, where a superconducting instability 
takes over at low temperatures (T < T c m 1.2 K). This superconducting state is very 
likely to be triplet, as sug gest ed by upper critical fields and NMR measurements. 
Recent experiments 0, liol lilt suggest a region of macroscopic coexistence between TSC 
and SDW, where both orders are non-uniform. This coexistence region can be related 
to existing theoretical proposals. For instance, strictly one-dimensional theories invoking 
SO (4) symmetry 2], or negative interface energies have allowed for coexisting TSC 
and SDW. However, these previous theories are not directly applicable to three-dimensional 
but highly anisotropic superconductors like the Bechgaard salts, where the SO (4) symmetry 
is absent, and negative interface energies are not necessary conditions for the coexistence. In 
this manuscript, we derive microscopically the pressure versus temperature phase diagram 
indicating the TSC, the SDW and the TSC/SDW phases, and show that the TSC and 
SDW order parameters are both non-uniform in the coexistence region. Furthermore, we 
investigate the effects of an external magnetic field and suggest that a canting transition 
of the SDW order parameter may occur, and alter the nature of the TSC state in the 
coexistence region. 
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The compound (TMTSF^PFg can be described approximately by an orthorhombic lattice 
with dispersion 

£k = ~\tx\ cos(k x a) — \t y \ cos(kyb) — \t z \ cos(k z c), (1) 

where transfer integrals \t x \, \t y \ and \t z \ satisfy the relations ^> \t y \ ^> \t z \ representing 
the quasi-one- dimensionality. Here, a, b, and c correspond to unit cell lengths along the 
crystallographic axes a(x), b'(y) and c*(z), respectively 

We use natural units (h = ks = c = 1) and work with Hamiltonian 7i = l Ho + 'H\- a t ) where 
the non-interacting part is Ho = a £k c iL a c k,cn an d £k = e k — A* is the dispersion shifted 
by the chemical potential, which may include a Hartree shift. The interaction part is 

n** = EE v ^ k ') d U k > p) ■ d ^( k '' p) 

kk'p aft-yS 

+ EE ^(q)s^(k,q)-s 7 ,(k',q), (2) 

kk'q a/37<5 

where the first and second terms describe interactions in TSC and SDW channels, respec- 
tively. These interactions allow for the possibility of competition or coexistence of TSC and 
SDW instabilities at low temperatures. Here, a, /3, 7 and 6 are spin indices and k, k', p 
and q represent linear momenta. The vector operator d^k, p) = c] c+ , 2 Q, v a/3 c ik+p/2 /3> anc ^ 
s^(k, q) = Cj t _ q / 2 a cr Q , / 3Ck+ q /2,/3- The matrix v a/ 3 = {itTa y ) a p, and <Tj are Pauli matrices. In 
the case of weak spin-orbit coupling, the TSC interaction V(k, k') can be choosen as 

V(k, k') = Vh r (k, k')0r(k)0r(k') (3) 

where V is a prefactor with dimension of energy, /ir(k, k') [</>p(k)] characterizes the mo- 
mentum dependence [symmetry basis function] for an irreducible representation T of the 
orthorhombic group D 2 h- Without loss of the generality regarding symmetry properties, we 
take /ir(k, k') = 1, and consider only unitary triplet states corresponding to p x symmetry. 
The order parameter for TSC can be defined as 

D(p) = (^^ r (k)d a/3 (k,p)), (4) 

k,a/3 

while the SDW order parameter can be defined as 

N(q) = J(q)<5> Q/3 (k, q )>. (5) 

M/3 
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With these definitions, the effective Hamiltonian is 

Hcs = Hq + Htsc + T^sdw, (6) 

where the TSC contribution is Htsc = EpMpVEk,^ 0r(k)d a/3 (k, p)+H.C.]-£ p Dt(p). 
D(p)/y, and the TSC term is H S dw = E q t N (-q) ' E k , Q/3 S a/3 (k, q) + H.C.] - £ q N(-q) • 
N(q)/J(q). The effective action of this Hamiltonian as a function of D(p) and N(q) is 
obtained by integrating out the fermions. The quadratic terms are 

S 2 TSC = 5>(p)Dt(p).D(p), 
p 

5JDW = ^ S (q)N(-q).N(q), (7) 
q 

with coefficients A(p) = - 2T^ fc G(k + p/2,c; n )0r(k + p/2)G(-k + 

p/2,-u; n )0r(-k + p/2), an d B(q) = -2 J-^q) - 2T £ fe G(k, uj n )G(k + q, w ft ). Here, T 
is the temperature, G(k, u) n ) = l/(iu n + £k) is the fermion propagator, oj n is fermionic 
Matsubara frequency, and J2k = Skw„ ^ s use d to shorthand notation. 

In what follows we make some standard assumptions. First, we assume that the saddle 
point TSC order parameter is dominated by the zero center of mass momentum component 
D = D(p = 0). Second, we assume that the saddle point SDW order parameter N is a real 
vector in r-space, and that it has Fourier components determined by Fermi surface nesting 
vectors q = Qj = (±Q a , ±Qb, ±Q c ) 

. In this case, the coefficients -B(Qj) are identical for 
all Qj's, since the lattice dispersion in invariant under reflections and inversions compatible 
with the _D2h group. In addition, the coefficients of all higher order terms involving N(Qj)'s 
share the same properties. Given that N(r) is real, and that we have periodic boundary 
conditions, we can choose a specific reference phase where N(Qj) are real and identical. 
Thus, we define N = N(Q$) for all i, and the quadratic terms are dominated in the long 
wavelength limit by S^ sc w A(0)|D | 2 and Sf DW w (m/2) J B(Qi)|N | 2 , respectively. Here, 
m is the number of nesting vectors, and Q x = (Q a , Q b , Q c ) is chosen for definiteness. 

Notice that the two order parameters D(p) and N(q) do not couple to quadratic or- 
der, because TSC and SDW are instabilities in particle-particle and particle-hole channels, 
respectively. Thus, the two orders are independent to this order, and their corresponding 
vector order parameters are free to rotate. However, this freedom is lost when fourth order 
terms are included. 
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The coupling between D and N in fourth order is given by 

Sf = (Ci + C 2 /2)|D | 2 |N | 2 - C 2 \B • N | 2 , (8) 

where G\ = mT J2k u n )G 2 (— k, — u n )G(— k + Q 1; — cj n )0 r (— k)0 r (k), and the other co- 
efficient C 2 = mTJ2kG(k,oj n )G(\i+Q ll u n )G(—'k-Qi,uj n )G(—'k, —uj n )<j) r (— k — Q x )0 r (k). 
Notice that the second term on Eq. (JHJ) can be parametrized as — C 2 cos 2 (#)|Do| 2 |No| 2 , where 
cos 2 6 = |D ■ N | 2 /|D | 2 |N | 2 < 1 is independent of |D | and |N |. Since D is unitary, its 
global phase can be eliminated in S±, and 9 can be regarded as the angle between D and 
No. The coefficient C 2 for (TMTSF) 2 PF 6 is positive, indicating that D and N tend to be 
aligned {6 = 0) or anti-aligned (9 — ir). 
Additional fourth order terms are 

Sj sc = I> a |Doh St™ = D 2 \N \\ (9) 

where D 1 = T J2 k G 2 (k, u n )G 2 (-k, -cu n )0 2 (k)0 2 (-k), and D 2 
(m/2)T E fc G(k, ^„)G 2 (k + Q 1; co n ) [G(k, u n ) + G(k+2Q 1 , u n )} + (m/4)T J2 k G(k, ^)G(k + 
Qi,^n) Z)i G ( k + Qi + Qi,^n)[C(k + Qj,o; n ) + G(k + Q 1 ,u n )}. Here Yli represents 
Y2q ^±Qj- This leads to the effective action 

S eS = S + S 2 + S 4 , (10) 

where Sq is the normal state contribution, S 2 = £j sc + Sf DW , and S4 = Di|D | 4 + -D2|N | 4 -|- 
C(6 ) )|D | 2 |N | 2 with C{6) = C x + C 2 /2 — C 2 cos 2 #. The phase diagram that emerges from 
this action leads to either bicritical or tetracritical points as illustrated in Fig. ^ When 
R = C 2 (0)/(ADiD 2 ) > 1 the critical point (P C ,T C ) is bicritical and there is a first order 
transition line at (m/2)5(Qi) = A(0) when both -B(Qi) < and A(0) < 0, as seen in 
Fig. da). However when R < 1, (P C ,T C ) is tetracritical and a coexistence region for 
TSC and SDW orders occurs when both £>(Qi) < and A(0) < 0, as seen in Fig. HJ^b). 
The action S e g obtained in three dimensions is not SO(4) invariant, and SO(4) symmetry 
based theories |12] can only be applied to one-dimensional systems, but not to the highly 
anisotropic but three-dimensional Bechgaard salts. 

The ratio R w (0.12) 2 < 1 for the Bechgaard salt (TMTSF) 2 PF 6 around (P C ,T C ), when 
the interaction strengths V, J are chosen to give the same T c = 1.2 K at quarter filling for 
parameters \t x \ = 5800 K, \t y \ = 1226 K, \t z \ = 58 K, used in combination with 0r(k) = 
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FIG. 1: P-T phase diagrams indicating (a) first order transition line with no coexistence phase, 
and (b) two second order lines with coexistence region between TSC and SDW phases. 

sm(k x a) (p^-symmetry for TSC) and the nesting vectors Q = (±7r/2a, ±7r/26, 0) (m = 
4). This shows that (TMTSF^PFr has an TSC/SDW coexistence region as suggested by 



experiments |9 
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11] 



To investigate the TSC/SDW coexistence region the effective action (JTU|) [with Q = 
(±Q a , 0)] is Fourier transformed into real space to give the Ginzburg-Laudau (GL) 
free energy density 

F — 3~n + -^TSC + -^SDW + ^C) (H) 

where T n is the normal state contribution, and Tq = C(0)|N(r)| 2 |D(r)| 2 is the coupling term 
of the two order parameters. The TSC contribution is .Ftsc = y 4(0)|D(r)| 2 + -Di|D(r)| 4 + 
Ylij 7Tsc[^-^( r )] ' [^i-D(r)], where 7tsc i s obtained from a small p expansion of ^4(p). The 
SDW contribution is J^sdw = / dr' [B(r, r')N(r) ■ N(r')] + (D 2 /m 2 )\N(v)\\ where B(r, r') is 
the Fourier transform of B(q) in Eq. (J7J). For the Bechgaard salt parameters, the prefactor 
C(0) of the coupling term Tq is positive, and hence represents a local repulsive interaction 
between the TSC and SDW order parameters. As a consequence, the TSC order parameter 
is non-uniform in the TSC/SDW coexistence region, and has a modulation induced by the 
SDW order parameter. Since R <C 1 for (TMTSF) 2 PF 6 , the coupling term Tq is small 
in comparison with the other fourth order coefficients Di,D 2 , and a perturbative solution 
is possible for |D(r)| and |N(r)|. At assumed zero TSC/SDW coupling C(0) = 0, the 
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FIG. 2: Magnitude square of (a) TSC and (b) SDW order parameters in the coexistence region. 



saddle point modulation for the SDW order parameter is N(r) = mN cos(Qi • r), with 
| No | = [— m5(Qi)/3Z^2] 1 ^ 2 , while the saddle point for the TSC order parameter is D(r) = D 
with |D | = [-A(0)/2D 1 ] 1 / 2 . Including the coupling T G the new solution for the magnitude 
of TSC order parameter is 



liVr 
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|D(r)|-|D | = -VTjrrR 1/2 



cos(2Q a x) . cos(2Q b y) 



+ 



cos(2Q a a:) cos(2Q b y) 1 
+ 4 + + 8QQ 2 b 4 



(12) 



which shows explicitly 2Q a and 2Q b modulations along the a and b' axes, respectively. 



Here, & = \/|7tscM(0)I represents the TSC coherence length along the % direction, and 
v = (6D 2 /D 1 ) 1 / 2 . Notice that the modulation in |D(r)| disappears as the SDW order 
goes away |N | — > 0. The qualitative behavior of |D(r)| is shown in Fig. |2Ja). The new 
solution for the SDW order parameter to the first order correction is |N(r)| = — 
i? 1 / 2 |D | 2 /4v|No| 2 )|N | cos(Qj • r), and can be seen in Fig. Efb). Notice that the maxima 
of |D(r)| coincide with the minima of |N(r)| indicating that the TSC and SDW orders try 
to be locally excluded. Since the TSC and SDW modulations are out of phase, experiments 
that are sensitive to the spatial distribution of the spin density or Cooper pair charge density 
may reveal the coexistence of these inhomogeneous phases. 

Next, we analyze the effect of magnetic fields on this coexistence region. A uniform 
magnetic field H couples with charge via the Peierls substitution k — > k — |e|A in the 
dispersion relation given in Eq. (0), where A is the vector potential, and couples with spin 
via the paramagnetic term Hp = — //qH • a/3 c^a-a^c^, where /i is the effective magnetic 



7 



moment. Thus, the effective Hamiltonian becomes 



e| A) + Htsc + ^sdw + Hp- 



(13) 



Upon integration of the fermions, the corresponding effective action is 



5 eff (H) = 5 (H) + 5 2 (H) + 5 4 (H), 



(14) 



where 5q(H) = 5*0 + |H| 2 /87r — x„|H| 2 /2, Xn is the uniform electronic spin suscepti- 
bility of the normal state, 5 2 (H) is obtained from 5*2 by the Peierls substitution, and 
S4(H) = 5 4 + {Eh. + £ 2 /2)|H| 2 |D | 2 - E 2 \H ■ D | 2 + (Ft - F 2 /2)|H| 2 |N | 2 + F 2 |H ■ 
N | 2 . The coefficients are E x = 2/i 2 T £ fc G 3 (-k, uj n )G(k, -w n )0 r (k)0r(-k), E 2 = 
2fi 2 TJ2 k G 2 (k,uj n )G\-k,-uj n )Mk)M-k), Fx = m/i 2 T£ fc G 3 (k,^ n )[G(k + Q u u n ) + 
G(k — Qi, u> n )], and F 2 = m/i^T J2k ^ 2 (k, ^n)C 2 (k + Qi, u n ). A detailed calculation shows 
that the coefficient E\ = —E 2 /2, hence the coupling of H to D can be described in the more 
familiar form F M - where = Xn^ + E 2 D*^D U . 

For Bechgaard salts, the coefficients E 2 < and F 2 > indicating that D and N prefer 
to be perpendicular to the magnetic field H. These conditions, when combined with C 2 > 
in Eq. (jHJ), indicate that D and N prefer to be parallel to each other, but perpendicular 
to H. However, the relative orientation of these vectors in small fields is affected by spin 
anisotropy effects which were already observed in (TMTSF) 2 PF 6 , where the easy axis for N 
is the b' direction [(|. Such an anisotropy effect can be described by adding a quadratic term 
—unN 2 , with un > 0, which favors N || b'. Similarly, the D vector also has anisotropic effect 
caused by spin-orbit coupling, and can be described by adding a quadratic term —u^Df, 
where % is the easy axis for TSC. (Quartic TSC and SDW terms also become anisotropic.) 

However, a sufficiently large H || b' can overcome spin anisotropy effects, and drive the 
N vector to flop onto the a-c* plane. This canting (flop) transition was reported (f| in 
(TMTSF) 2 PF 6 for H « 1 T at zero pressure and T = 8 K. If such a spin-flop transition 
persists near the TSC/SDW critical point (P c , T c ) as suggested in our discussion, then the 
flop transition of the N vector forces the D vector to flop as well, and has potentially serious 
consequences to the superconducting state. Schematic phase diagrams are are shown in 
Fig. 3(a) and 3(b). For P < P c , if a flop transition occurs for Hp < -Hi(O) [see Fig. 3(a)], 
then N flops both in the pure SDW and in the TSC/SDW coexistence phases, in which case it 
forces D vector to flop as well. If the flop transition occurs for i7g DW (0) < H F < i?i(0) (not 
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FIG. 3: H-T phase diagrams showing the TSC/SDW coexistence region (thick solid line) and 
canting transitions (double line) for (a) P < P c and (b) P > P c . 

shown) then only the pure SDW phase is affected. This situation is qualitatively different 
for P > P c . In the zero (weak) spin-orbit coupling limit the D vector is free to rotate in a 
magnetic field and tends to be perpendicular to H in order to minimize its magnetic free 
energy Fm- Thus, for H || b' and |H| > H 2 , the D vector lies in the a-c* plane since there 
is no SDW order. However, at lower temperatures and small magnetic fields when TSC and 
SDW orders coexist, the spin anisotropy field forces N to be along b r / and N forces D to 
flop from the a-c* plane to b' direction. This canting transition occurs at Hp < #2(0) [see 
Fig. 3(b)], when N flops in the TSC/SDW coexistence phase, and forces the D vector to 
flop as well. 

In summary, we showed that the TSC and SDW order parameters can coexist in the 
P-T phase diagram of quasi-one- dimensional organic conductors. In the coexistence region 
the TSC order parameter is non-uniform, and its modulation is induced via the SDW order 
parameter. We also showed that theories based on SO (4) symmetry cannot be applied to 
these highly anisotropic three-dimensional systems, since they are strictly valid only in the 
one-dimensional limit. Furthermore, we discussed qualitatively magnetic field effects on the 
coexistence region. We proposed that a magnetic field induced canting transition of the SDW 
order parameter affects dramatically the phase diagram of the coexistence region, both below 
and above the critical pressure. We would like to thank NSF for support (DMR-0304380). 
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